Pangenome graphs improve the analysis of structural variants in rare genetic diseases.

Rare DNA alterations that cause heritable diseases are only partially resolvable by clinical next-generation sequencing due to the difficulty of detecting structural variation (SV) in all genomic contexts. Long-read, high fidelity genome sequencing (HiFi-GS) detects SVs with increased sensitivity and enables assembling personal and graph genomes. We leverage standard reference genomes, public assemblies (n = 94) and a large collection of HiFi-GS data from a rare disease program (Genomic Answers for Kids, GA4K, n = 574 assemblies) to build a graph genome representing a unified SV callset in GA4K, identify common variation and prioritize SVs that are more likely to cause genetic disease (MAF < 0.01). Using graphs, we obtain a higher level of reproducibility than the standard reference approach. We observe over 200,000 SV alleles unique to GA4K, including nearly 1000 rare variants that impact coding sequence. With improved specificity for rare SVs, we isolate 30 candidate SVs in phenotypically prioritized genes, including known disease SVs. We isolate a novel diagnostic SV in KMT2E, demonstrating use of personal assemblies coupled with pangenome graphs for rare disease genomics. The community may interrogate our pangenome with additional assemblies to discover new SVs within the allele frequency spectrum relevant to genetic diseases.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

Single duplex DNA sequencing with CODEC detects mutations with high sensitivity.

Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next-generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules ('single duplexes') to discern the true mutations on both strands. Here we present Concatenating Original Duplex for Error Correction (CODEC), which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than duplex sequencing. CODEC revealed mutation frequencies of 2.72 × 10-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide, clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability with 10-fold greater sensitivity and mutational signatures, and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations, which are commonly obscured by NGS errors.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing.

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Biochemical characterization of two novel mutations in the human high-affinity choline transporter 1 identified in a patient with congenital myasthenic syndrome.

Congenital myasthenic syndrome (CMS) is a heterogeneous condition associated with 34 different genes, including SLC5A7, which encodes the high-affinity choline transporter 1 (CHT1). CHT1 is expressed in presynaptic neurons of the neuromuscular junction where it uses the inward sodium gradient to reuptake choline. Biallelic CHT1 mutations often lead to neonatal lethality, and less commonly to non-lethal motor weakness and developmental delays. Here, we report detailed biochemical characterization of two novel mutations in CHT1, p.I294T and p.D349N, which we identified in an 11-year-old patient with a history of neonatal respiratory distress, and subsequent hypotonia and global developmental delay. Heterologous expression of each CHT1 mutant in human embryonic kidney cells showed two different mechanisms of reduced protein function. The p.I294T CHT1 mutant transporter function was detectable, but its abundance and half-life were significantly reduced. In contrast, the p.D349N CHT1 mutant was abundantly expressed at the cell membrane, but transporter function was absent. The residual function of the p.I294T CHT1 mutant may explain the non-lethal form of CMS in this patient, and the divergent mechanisms of reduced CHT1 function that we identified may guide future functional studies of the CHT1 myasthenic syndrome. Based on these in vitro studies that provided a diagnosis, treatment with cholinesterase inhibitor together with physical and occupational therapy significantly improved the patient's strength and quality of life.

Serial enrichment of heteroduplex DNA using a MutS-magnetic bead system.

Numerous applications in molecular biology and genomics require characterization of mutant DNA molecules present at low levels within a larger sample of non-mutant DNA. This is often achieved either by selectively amplifying mutant DNA, or by sequencing all the DNA followed by computational identification of the mutant DNA. However, selective amplification is challenging for insertions and deletions (indels). Additionally, sequencing all the DNA in a sample may not be cost effective when only the presence of a mutation needs to be ascertained rather than its allelic fraction. The MutS protein evolved to detect DNA heteroduplexes in which the two DNA strands are mismatched. Prior methods have utilized MutS to enrich mutant DNA by hybridizing mutant to non-mutant DNA to create heteroduplexes. However, the purity of heteroduplex DNA these methods achieve is limited because they can only feasibly perform one or two enrichment cycles. We developed a MutS-magnetic bead system that enables rapid serial enrichment cycles. With six cycles, we achieve complete purification of heteroduplex indel DNA originally present at a 5% fraction and over 40-fold enrichment of heteroduplex DNA originally present at a 1% fraction. This system may enable novel approaches for enriching mutant DNA for targeted sequencing.

An HNRNPK-specific DNA methylation signature makes sense of missense variants and expands the phenotypic spectrum of Au-Kline syndrome.

Au-Kline syndrome (AKS) is a neurodevelopmental disorder associated with multiple malformations and a characteristic facial gestalt. The first individuals ascertained carried de novo loss-of-function (LoF) variants in HNRNPK. Here, we report 32 individuals with AKS (26 previously unpublished), including 13 with de novo missense variants. We propose new clinical diagnostic criteria for AKS that differentiate it from the clinically overlapping Kabuki syndrome and describe a significant phenotypic expansion to include individuals with missense variants who present with subtle facial features and few or no malformations. Many gene-specific DNA methylation (DNAm) signatures have been identified for neurodevelopmental syndromes. Because HNRNPK has roles in chromatin and epigenetic regulation, we hypothesized that pathogenic variants in HNRNPK may be associated with a specific DNAm signature. Here, we report a unique DNAm signature for AKS due to LoF HNRNPK variants, distinct from controls and Kabuki syndrome. This DNAm signature is also identified in some individuals with de novo HNRNPK missense variants, confirming their pathogenicity and the phenotypic expansion of AKS to include more subtle phenotypes. Furthermore, we report that some individuals with missense variants have an "intermediate" DNAm signature that parallels their milder clinical presentation, suggesting the presence of an epi-genotype phenotype correlation. In summary, the AKS DNAm signature may help elucidate the underlying pathophysiology of AKS. This DNAm signature also effectively supported clinical syndrome delineation and is a valuable aid for variant interpretation in individuals where a clinical diagnosis of AKS is unclear, particularly for mild presentations.

Applications of Single-Cell DNA Sequencing.

Over the past decade, genomic analyses of single cells-the fundamental units of life-have become possible. Single-cell DNA sequencing has shed light on biological questions that were previously inaccessible across diverse fields of research, including somatic mutagenesis, organismal development, genome function, and microbiology. Single-cell DNA sequencing also promises significant future biomedical and clinical impact, spanning oncology, fertility, and beyond. While single-cell approaches that profile RNA and protein have greatly expanded our understanding of cellular diversity, many fundamental questions in biology and important biomedical applications require analysis of the DNA of single cells. Here, we review the applications and biological questions for which single-cell DNA sequencing is uniquely suited or required. We include a discussion of the fields that will be impacted by single-cell DNA sequencing as the technology continues to advance.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features.

PURPOSE: The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10. METHODS: Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry. RESULTS: A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein. CONCLUSION: We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

Complex Autoinflammatory Syndrome Unveils Fundamental Principles of JAK1 Kinase Transcriptional and Biochemical Function.

Autoinflammatory disease can result from monogenic errors of immunity. We describe a patient with early-onset multi-organ immune dysregulation resulting from a mosaic, gain-of-function mutation (S703I) in JAK1, encoding a kinase essential for signaling downstream of >25 cytokines. By custom single-cell RNA sequencing, we examine mosaicism with single-cell resolution. We find that JAK1 transcription was predominantly restricted to a single allele across different cells, introducing the concept of a mutational "transcriptotype" that differs from the genotype. Functionally, the mutation increases JAK1 activity and transactivates partnering JAKs, independent of its catalytic domain. S703I JAK1 is not only hypermorphic for cytokine signaling but also neomorphic, as it enables signaling cascades not canonically mediated by JAK1. Given these results, the patient was treated with tofacitinib, a JAK inhibitor, leading to the rapid resolution of clinical disease. These findings offer a platform for personalized medicine with the concurrent discovery of fundamental biological principles.

Transcriptome sequencing identifies a noncoding, deep intronic variant in CLCN7 causing autosomal recessive osteopetrosis.

BACKGROUND: Over half of children with rare genetic diseases remain undiagnosed despite maximal clinical evaluation and DNA-based genetic testing. As part of an Undiagnosed Diseases Program applying transcriptome (RNA) sequencing to identify the causes of these unsolved cases, we studied a child with severe infantile osteopetrosis leading to cranial nerve palsies, bone deformities, and bone marrow failure, for whom whole-genome sequencing was nondiagnostic. METHODS: We performed transcriptome (RNA) sequencing of whole blood followed by analysis of aberrant transcript isoforms and osteoclast functional studies. RESULTS: We identified a pathogenic deep intronic variant in CLCN7 creating an unexpected, frameshifting pseudoexon causing complete loss of function. Functional studies, including osteoclastogenesis and bone resorption assays, confirmed normal osteoclast differentiation but loss of osteoclast function. CONCLUSION: This is the first report of a pathogenic deep intronic variant in CLCN7, and our approach provides a model for systematic identification of noncoding variants causing osteopetrosis-a disease for which molecular-genetic diagnosis can be pivotal for potentially curative hematopoietic stem cell transplantation. Our work illustrates that cryptic splice variants may elude DNA-only sequencing and supports broad first-line use of transcriptome sequencing for children with undiagnosed diseases.

Correction: DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

Following the publication of the article, it was noted that the last column in Table 1, the total % should have read 5/8 (62.5) for the 'Epilepsy' row, and not 5.7 (71.4). This has now been amended in the HTML and PDF of the original article.

DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

DPH1 variants have been associated with an ultra-rare and severe neurodevelopmental disorder, mainly characterized by variable developmental delay, short stature, dysmorphic features, and sparse hair. We have identified four new patients (from two different families) carrying novel variants in DPH1, enriching the clinical delineation of the DPH1 syndrome. Using a diphtheria toxin ADP-ribosylation assay, we have analyzed the activity of seven identified variants and demonstrated compromised function for five of them [p.(Leu234Pro); p.(Ala411Argfs*91); p.(Leu164Pro); p.(Leu125Pro); and p.(Tyr112Cys)]. We have built a homology model of the human DPH1-DPH2 heterodimer and have performed molecular dynamics simulations to study the effect of these variants on the catalytic sites as well as on the interactions between subunits of the heterodimer. The results show correlation between loss of activity, reduced size of the opening to the catalytic site, and changes in the size of the catalytic site with clinical severity. This is the first report of functional tests of DPH1 variants associated with the DPH1 syndrome. We demonstrate that the in vitro assay for DPH1 protein activity, together with structural modeling, are useful tools for assessing the effect of the variants on DPH1 function and may be used for predicting patient outcomes and prognoses.

Integrated genome and transcriptome sequencing identifies a noncoding mutation in the genome replication factor DONSON as the cause of microcephaly-micromelia syndrome.

While next-generation sequencing has accelerated the discovery of human disease genes, progress has been largely limited to the "low hanging fruit" of mutations with obvious exonic coding or canonical splice site impact. In contrast, the lack of high-throughput, unbiased approaches for functional assessment of most noncoding variants has bottlenecked gene discovery. We report the integration of transcriptome sequencing (RNA-seq), which surveys all mRNAs to reveal functional impacts of variants at the transcription level, into the gene discovery framework for a unique human disease, microcephaly-micromelia syndrome (MMS). MMS is an autosomal recessive condition described thus far in only a single First Nations population and causes intrauterine growth restriction, severe microcephaly, craniofacial anomalies, skeletal dysplasia, and neonatal lethality. Linkage analysis of affected families, including a very large pedigree, identified a single locus on Chromosome 21 linked to the disease (LOD > 9). Comprehensive genome sequencing did not reveal any pathogenic coding or canonical splicing mutations within the linkage region but identified several nonconserved noncoding variants. RNA-seq analysis detected aberrant splicing in DONSON due to one of these noncoding variants, showing a causative role for DONSON disruption in MMS. We show that DONSON is expressed in progenitor cells of embryonic human brain and other proliferating tissues, is co-expressed with components of the DNA replication machinery, and that Donson is essential for early embryonic development in mice as well, suggesting an essential conserved role for DONSON in the cell cycle. Our results demonstrate the utility of integrating transcriptomics into the study of human genetic disease when DNA sequencing alone is not sufficient to reveal the underlying pathogenic mutation.

Resolving rates of mutation in the brain using single-neuron genomics.

Whether somatic mutations contribute functional diversity to brain cells is a long-standing question. Single-neuron genomics enables direct measurement of somatic mutation rates in human brain and promises to answer this question. A recent study (Upton et al., 2015) reported high rates of somatic LINE-1 element (L1) retrotransposition in the hippocampus and cerebral cortex that would have major implications for normal brain function, and suggested that these events preferentially impact genes important for neuronal function. We identify aspects of the single-cell sequencing approach, bioinformatic analysis, and validation methods that led to thousands of artifacts being interpreted as somatic mutation events. Our reanalysis supports a mutation frequency of approximately 0.2 events per cell, which is about fifty-fold lower than reported, confirming that L1 elements mobilize in some human neurons but indicating that L1 mosaicism is not ubiquitous. Through consideration of the challenges identified, we provide a foundation and framework for designing single-cell genomics studies.

Somatic mutation in single human neurons tracks developmental and transcriptional history.

Neurons live for decades in a postmitotic state, their genomes susceptible to DNA damage. Here we survey the landscape of somatic single-nucleotide variants (SNVs) in the human brain. We identified thousands of somatic SNVs by single-cell sequencing of 36 neurons from the cerebral cortex of three normal individuals. Unlike germline and cancer SNVs, which are often caused by errors in DNA replication, neuronal mutations appear to reflect damage during active transcription. Somatic mutations create nested lineage trees, allowing them to be dated relative to developmental landmarks and revealing a polyclonal architecture of the human cerebral cortex. Thus, somatic mutations in the brain represent a durable and ongoing record of neuronal life history, from development through postmitotic function.

Loss of PCLO function underlies pontocerebellar hypoplasia type III.

OBJECTIVE: To identify the genetic cause of pontocerebellar hypoplasia type III (PCH3). METHODS: We studied the original reported pedigree of PCH3 and performed genetic analysis including genome-wide single nucleotide polymorphism genotyping, linkage analysis, whole-exome sequencing, and Sanger sequencing. Human fetal brain RNA sequencing data were then analyzed for the identified candidate gene. RESULTS: The affected individuals presented with severe global developmental delay and seizures starting in the first year of life. Brain MRI of an affected individual showed diffuse atrophy of the cerebrum, cerebellum, and brainstem. Genome-wide single nucleotide polymorphism analysis confirmed the linkage to chromosome 7q we previously reported, and showed no other genomic areas of linkage. Whole-exome sequencing of 2 affected individuals identified a shared homozygous, nonsense variant in the PCLO (piccolo) gene. This variant segregated with the disease phenotype in the pedigree was rare in the population and was predicted to eliminate the PDZ and C2 domains in the C-terminus of the protein. RNA sequencing data of human fetal brain showed that PCLO was moderately expressed in the developing cerebral cortex. CONCLUSIONS: Here, we show that a homozygous, nonsense PCLO mutation underlies the autosomal recessive neurodegenerative disorder, PCH3. PCLO is a component of the presynaptic cytoskeletal matrix, and is thought to be involved in regulation of presynaptic proteins and synaptic vesicles. Our findings suggest that PCLO is crucial for the development and survival of a wide range of neuronal types in the human brain.